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We study the roles of the dynamical high order perturbation and statistically non-linear infrared 
fluctuation/correlation in the virial equation of state for the Fermi gas in the unitary limit. In- 
corporating the quantum level crossing rearrangement effects, the spontaneously generated entropy 
departing from the mean-field theory formalism leads to concise thermodynamical expressions. The 
dimensionless virial coefficients with complex non-local correlations are calculated up to the fourth 
order for the first time. The virial coefficients of unitary Fermi gas are found to be proportional to 
those of the ideal quantum gas with integer ratios through a general term formula. Counterintu- 
itively, contrary to those of the ideal bosons (a^ = — ) or fermions^ ' = ^75), the second virial 
coefficient a-2 of Fermi gas at unitarity is found to be equal to zero. With the vanishing leading order 
quantum correction, the BCS-BEC crossover thermodynamics manifests the famous pure classical 
Boyle's law in the Boltzmann regime. The non-Gaussian correlation phenomena can be validated 
by studying the Joule- Thomson effect. 
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I. INTRODUCTION 

Unconventional unitary fermion physics is associated 
with a variety of strongly interacting topics. This theme 
can test the many-body theories, from neutron stars and 
nuclear matter to quark-gluon plasmas etc. 

In recent years, the study for the interacting fermion 
matter properties has attracted much attention in quan- 
tum many-body community[l]. This is attributed to the 
rapid progress of atomic Fermi gas experiments. Control- 
ling the 5-wavc scattering length between two different 
spin components allows one to control the interaction 
strength by using a magnetically tuned Feshbach reso- 
nance. With the magnetic tuning technique, increasing 
the interaction strength of atomic fermions with scatter- 
ing length a from — oo to +oo resulting in bound bo- 
son systems exhibits the Bardeen-Cooper-Schrieffer to 
Bose-Einstein condensation (BCS-BEC crossover). The 
two regimes with positive and negative S-wave scatter- 
ing length meet in the strongly interacting limit with 
the divergent scattering length. At the resonant point, 
the scattering cross-section will be saturated as a — 
4ir/k 2 (with k being the relative wave- vector magnitude 
of the colliding particles) due to the fundamental uni- 
tary property limit. The strongly interacting BCS-BEC 
crossover topic is literally called the unitary Fermi gas 
thermodynamics^— 16]. 

Generally, the thermodynamical properties of the low 
energy dilute Fermi system are determined by the S- 
wave scattering length a, the particle number density n 
as well as the temperature T. In the resonant regime 
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with a zero-energy bound state, the divergent scatter- 
ing length will certainly drop out in the thermodynam- 
ical quantities; i.e., the thermodynamical properties are 
universal [1-3]. The divergent scattering length poses an 
intractable many-body problem. 

In addressing the BCS-BEC crossover thermodynam- 
ics, the fundamental issue is the ground state energy. On 
the basis of the dimensional analysis, the dimensionless 
coefficient £ relates the energy per particle E/N = £§e/ 
with the Fermi kinetic energy e/ = fc?/(2m). Here, 
to is the bare fermion mass while kf is the Fermi mo- 
mentum. The fundamental universal coefficient £ has 
prompted many theoretical or experimental efforts in re- 
cent years[l]. Furthermore, the finite temperature ther- 
modynamic properties of unitary Fermi gas are as in- 
triguing as the zero-temperature ground state energy 
and many experimental/theoretical efforts have been 
made[5, 6] [7, 8]. 

Compared with the zero-temperature ground state en- 
ergy, the additional energy scale, i.e, the reciprocal ther- 
modynamical de Broglie wavelength A -1 = \J mT/(2n) 
complicates the theme and makes the universal prop- 
erty analysis more profound [3]. In the weak degenerate 
Boltzmann regime, the universal properties are charac- 
terized by the virial coefficients [3, 4, 9, 15]. For exam- 
ple, the virial equation of state is related with the neu- 
trinosphere physics in supernovae for the dilute nuclear 
matter matter. It is believed that the virial equation of 
state will influence the detailed information of the neu- 
trino response of low-density neutron matter [9]. How to 
calculate the virial coefficient in the unitary limit regime 
remains an important many-body topic. Like in address- 
ing the zero-temperature ground state energy problem, 
the central task is to understand the novel non-linear 
fluctuation/correlation physics. 
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The virial expansion is the basic tool for use in dis- 
cussing the thcrmodynamical properties and should be 
model independent. In thermodynamics, the non-linear 
virial expansion is the infinite scries of the pressure ac- 
cording to the particle number density. Even as a funda- 
mental theme, this question is very challenging and by no 
means resolved yet. To derive the high order virial coeffi- 
cients of strongly correlated fermions in the unitary limit, 
the involved quantum statistical fluctuation/correlation 
and detailed dynamical effects must be clarified, which 
is the novel systematic requirement for the sound theo- 
retical efforts. With the aim of calibrating the universal 
virial coefficients, we strive to examine the spontaneously 
quantum level shift contribution on the entropy (counting 
the microscopic states) in a dynamically and thermody- 
namically self-consistent way. 

In this work, the dynamical and statistical correlation 
analysis explicitly demonstrates that the dense and hot 
thermodynamics at unitarity obey the virial theorem for 
the ideal non-interacting gas, i.e., P = 2/3E/V[3, 5]. 
Meanwhile, the calculated virial coefficients at unitarity 
are found to be proportional to those for the ideal Fermi 
gas with integer ratios through an universal general term 
formula. 

In the strongly interacting system, the dynam- 
ical effects compete with the non-linear fluctua- 
tion/correlations. The second virial coefficient is found to 
be vanishing due to the complicated correlations. From 
the viewpoint of bulk properties, is the dilute unitary 
Fermi gas behavior much more like the classical ideal 
Boyle gas in the weak degenerate Boltzmann regime? 

The present paper is organized as follows. In Sec. II, 
the statistical method with effective field theory formal- 
ism is presented. The approach is promoted by the phase 
separation-instability discussions for a compact environ- 
ment containing a competitive Coulomb frustration ele- 
ment. In this work, we will take into account the novel 
non-Gaussian fluctuation effects on the strongly inter- 
acting fermions thermodynamics. The thermodynamical 
consistency, virial theorem discussion and the high order 
virial coefficient calculations of the unitary Fermi sys- 
tem are presented in Sec. III. The numerical results are 
also presented in this section in order to allow compari- 
son with previous investigations. In Sec. IV, the scaling 
property of equation of state and the virial expansion in 
the unitary regime are further analyzed. Summarizing 
remarks are made in Sec.V. 

The calculations are performed in terms of the uni- 
versal four-fermions contact interaction formalism. The 
natural units ks = h = 1 are used throughout the paper. 



II. STATISTICAL DYNAMICS WITH 
NON-GAUSSIAN CORRELATIONS 

The strongly interacting matter offers a plausible per- 
spective in looking for the general statistical field theory 
methods. In the presence of a medium, the unitary topic 



becomes very challenging because the system is strongly 
correlated and has no small parameter applicable for any 
controlled perturbative calculations. 

Essentially, the strong correlation effects are highly 
non-linear and the pronounced turbulent features will ap- 
pear in the unitary system. The conventional mean-field 
theory or loop diagram ring and ladder resummation per- 
turbative techniques cannot be employed for the unitary 
fermions. Consequently, more theoretical attempts are 
urgently needed to understand the detailed dynamical 
role of the strongly interacting fermions thermodynam- 
ics. Although there are tremendous updating efforts, a 
soundly exact theory concerning the behavior at unitar- 
ity is still not available. 

With the goal of calculating the high order virial coeffi- 
cients, the procedure developed can allow the systematic 
rearrangements of the individual expansions while avoid- 
ing the theoretical double countings. The thermodynam- 
ical expressions obtained are in parallel with those with 
the linear bare contact interaction formalism. Therefore, 
the method developed will be referred to as the quasi- 
Gaussian/quasi-linear approximation in order to indicate 
the difference from and similarity to lowest order mean- 
field theory. 



A. In-medium effective action 

The in-mcdium behavior associated with the many- 
body characteristic is the key, while a non-perturbative 
approach is crucial. In many-body theory, an established 
fundamental perspective is that particles can change 
their spectrum properties in a dense and hot strongly cor- 
related medium. These changes will be reflected in the 
mass shifts and/or in the development of excited complex 
spectral property modifications. Considering the coun- 
teracting influences of the surrounding medium and spon- 
taneous single-particle spectrum modification due to the 
off-shell dispersive effects, a medium-scaling functional 
has been proposed[17] 

U*Jn,T] r 
+ ^ J d 3 xr a (x)r (x)Mx)Mx). (1) 

In Eq.(f), a, (3 =|,| represent the (hyperfine-)spin pro- 
jection Ising- variable. The effective Hamiltonian is the 
same as the original Bethe-Peierls zero-range contact in- 
teraction version, except that the bare coupling constant 
Uq = 4ira/m is substituted by an effective medium- 
scaling functional U*q[u, T]. With the bare potential 
Uq and corresponding vanishing fi ra in the vacuum limit 
n — > 0, the Hamiltonian H reduces to the original version 
possessing a global U(l ) or Z2 gauge symmetry. 

Due to the explicit medium dependence of U*^ in 
Eq.(l), we introduce the additional counterterm STi. oc 
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H r [n,T], which is enforced by the fundamental thermo- 
dynamical Hugenholtz-van Hove (HvH) theorem[18, 19]. 
Without loss of generality, we take care of the fully sym- 
metric scenario with [i ra = \i r . The complementary SH 
implies that the correlation effects on the single parti- 
cle energy spectrum are further taken into account as 
an effective single-body potential or spontaneously gen- 
erated binding energy in the spirit of density functional 
theory [20]. On average, the thermodynamic vacuum will 
have been "shifted" by fi r N non-perturbatively. 

As remarked in Rcf.[17, 21], the central ingredient of 
this non-perturbative procedure is the way of deriving 
the correlating coupling functional U*q. Essentially, its 
constitution is beyond the bare Bethe-Peierls dynamics 
itself; i.e, the derivation of U*q must be based on a more 
underlying physical law[18]. The functional U*g should 
reasonably encode the non-perturbative counteracting ef- 
fects of the surrounding environment. 

In thermodynamics, the surrounding environment 
plays a counteracting frustration role according to Le 
Chatelier's stability principle. This general principle ac- 
counts for the environment preventing an instantaneous 
departure from equilibrium with an alternating minus 
function[17, 22], which is consistent with the second law 
of thermodynamics. We find that the frustrating cor- 
relation effects of the surrounding environment can be 
realized via a twisted composite rearrangement matrix- 
vertex [17] 



U. 
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eff l-(m%/2)U ' 
I 



a cff 



1 — {2-Kw? D jm)a 



(2) 
(3) 



In order to making the analytical expressions as concise 
as possible, the inverse scattering length notation em- 
ployed, a e g, is defined in terms of U*g = Ana^/m. 

The reader will have noticed that there is an alter- 
nating "negative" sign difference in the denominator of 
Eq.(2) compared with the loop ring and ladder resum- 
mation perturbative techniques; i.e., the medium-scaling 
potential Eq.(2) appears as an instantaneously anti- 
screening formalism. This specific minus sign leads to 
quite different physical motivations and calculational de- 
tails. The non-linear screening formalism makes it possi- 
ble for us to incorporate the intermediate off-shell effects 
in an analytical way with 4-momcntum independent alge- 
bra equations, i.e., instead of numerically solving the var- 
ious coupled integral equations for the multi-points cor- 
relation Green functions. This refreshing attempt is mo- 
tivated by the particular conformal analogy and/or dis- 
cursion of the zero-range unitary Fermi problem with the 
universal instantaneous Coulomb correlation discussions 
in a compact confinement environment [23, 24] based on 
the relativistic continuum Dirac field theory formalism 
[25, 26], where involves the complex oscillatory instabili- 
ties with short range and long-range force's competitions. 



By analogy to the generalized Dyson- Schwinger calcu- 
lations with finite temperature Green function theory [17, 
21], but more conveniently, the Debye mass parameter 
m 2 D in Eq.(2) can be alternatively given by the general- 
ized Ward-Identity 



m 



D 



^L-TSJA/.M-W, (4 » 



where [i* is the effective chemical potential, the collec- 
tive implicit variable defined below. At T = 0, it reduces 
to the familiar m 2 D — kfm/ir 2 = 2N(ef), where N(ef) 
is the unperturbated density of states on the Fermi sur- 
face for one component fermions[27]. m 2 D characterizes 
the fluctuation physics because it is related to the well- 
known Pauli paramagnetic spin-spin or particle number 
susceptibility X — \ {dn/dfi) T (with an additional fac- 
tor | due to the degenerate degrees of freedom of the 
two-components symmetric system) according to 

(dn/d/i) T = (dn/dfi*) T (<9/x* /<9/i) T . (5) 



B. Grand thermodynamical potential with 
quasi-Gaussian approximation 

From the general Lagrange multiplier viewpoint, what 
we will perform is an evaluation of the relative min- 
imum fl(T,fl) of the shifted (0\H - {/J, - (i r )N\0) in- 
stead of directly evaluating the challenging absolute min- 
imum S}(T,/i) of the grand thermodynamic potential 
(0\H - juiV|0)[17, 21]. With auxiliary physical con- 
straints, the realistic f2(T, /x) incorporating the thermo- 
dynamical vacuum fluctuation and correlation effects is 
indirectly derived from the former for the given chemical 
potential \x. 

Firstly, the shifted relative minimum is evaluated 
by freezing the medium-dependent interaction potential 
with the conventional condensation formalism[17, 21] 



&(T, //*) 
V 



= —P + /j, r n 



na eS m 2 



n 2 -2T [ ln(l + e-' 3 (*^*)),(6) 

Jk 



with J, = Jc? 3 k/(27r) 3 being the momentum integral. 
The constraining self-consistent equation for the single 
particle Green function gives the definition of the effec- 
tive chemical potential /x* according to 



2na e« , * 
/i = /i — (i r = —n + n . 

m 



(7) 



The total number density n = nj + nj. = 2n| is given by 
2 // fc = n(T,/i*), (8) 

Jk 

with the defined quasi-particle Fermi-Dirac distribution 
functions 

h= ~ 4 , z' = e^. (9) 

z'-ie^ + 1 
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Different from the multiplier chemical potential ji, the 
dynamically collective variable \x* characterizes the ad- 
ditional correlation effects. With the chemical potential 
correction, the effective fugacity z' is analogous to the 
conventional one z = e^. Employing the quasi-particle 
Fermi-Dirac distribution functions Eq.(9), the thermody- 
namical formulae can manifest the standard Fermi inte- 
grals fj(z') with j = - \, \, ■ ■ ■ . 

Secondly, the remaining task is to determine the cor- 
rection term with physical constraints, from which the re- 
alistic grand thermodynamical potential will be uniquely 
determined. The shift terms oc fx r in the pressure and 
chemical potential are canceled out by each other in the 
Helmholtz free energy density / = F/V = —P + fin. 
However, the analytical expression for the relative shift 
strength /i r can be indirectly derived from the physical 
free energy density by relaxing the medium dependence 
of interaction potential. 

Invoking the thermodynamical relations, one can have 



P = n 



/i 



2 , d(fM 

dn 

dn 



(10) 
(11) 



Comparing the results derived from Eq.(10) and Eq.(ll) 
with Eq.(6) or Eq.(7), one has the definite analytical ex- 
pression of the relative shift strength 



Mr[n ,T]=C(T, M *)(^ff) n\ 
The employed correlation factor C is defined by 



C(T,fi*) =m D 



dm D \ I-\li\A 
dn ) T ~ 2Tf 1/2 [z>] 



(12) 



(13) 



D 



and C(T, /i*) are related to the high order den- 
sity/spin susceptibilities; i.e., they have the crystal clear 
physical connotations. With these two physical variables, 
the integrated grand partition function Eq.(6) or scaling 
equation of state and chemical potential Eq.(7) are re- 
duced to 



formalism via the frozen interaction; i.e., they have ex- 
actly the same structure of mean-field theory with linear- 
like interaction. The last curious shift term in Eq.(14a) 
non-Gaussianly proportional to the cubic particle num- 
ber density, oc n 3 , characterizes the high order contribu- 
tions beyond the former. This spontaneously generated 
term is picked up in a thermodynamical way by resuming 
the medium dependence of the interaction matrix. Phys- 
ically, this rearrangement effect considerably shifts the 
chemical potential of the particle distribution function 
as displayed by Eq.(14b). 

The correlation contributions manifested by the factor 
C arc explicitly combined with the dynamical high order 
ones. 

It is of crucial importance that the low and high order 
contributions are mixed with each other through a collec- 
tive variable, n*, effective chemical potential. One can 
see that the grand thermodynamical potential il(T,/j,) 
is not the naive polynomial expanded according to the 
bare vacuum interaction strength Uq. The dependence 
of the grand thermodynamical potential fl on the col- 
lective correlation variable fi* can be numerically elimi- 
nated in favor of the realistic physical chemical potential 
fi. It is physically pretty satisfactory that the equations 
Eq.(14a) and Eq.(14b) include the highly non-linear or 
turbulent correlation contributions. In the weak coupling 
limit, the equations readily reduce to those in terms of 
the lowest order mean-field theory, where the non-linear 
fluctuation/correlation contributions disappear. 



III. THERMODYNAMICS WITH QUANTUM 
REARRANGING CORRELATIONS 

We now examine the thermodynamical quantities with 
the self-consistent equations or the grand partition func- 
tion Eq.(14a) with Eq.(14b). The scenario we present 
is quite general although we have restricted ourselves to 
the low energy long-wavelength thermodynamics of the 
ideal fully symmetric system. For instance, the formal- 
ism and expressions can be easily extended to the asym- 
metric fermions phase separation analysis with unequal 
populations [28]. 
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with 



2T 



P ideal = ^3/5/2^') 



n d (14a) 
(14b) 

(15) 



which is similar to the ideal Fermi gas but with the ef- 
fective chemical potential as the collective variable. 

In Eq.(14a) and Eq.(14b), the first two terms appear 
as those obtained with the canonical Gaussian integral 



A. Thermodynamical consistencies, entropy and 
energy densities 

As discussed above, the non-Gaussian correlation ef- 
fects considerably complicate the one-one Legendre cor- 
responding relation of particle number N and the mul- 
tiplier chemical potential fj, in the strongly interacting 
systems. Due to the explicit high order corrections for 
the physical chemical potential, the calculations require 
additional great care; i.e., one must check if the thermo- 
dynamical relations are exactly ensured before making 
further discussions. 

For instance, an obvious check is to verify that the 
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partial derivative of P according to fj, gives again the 
particle number density n = N/V Eq.(8) 



dP 



dP 



dp,* 
d/j, 



^3/3/2(^.(16) 



Eq.(16) is one of the stringent thermodynamic con- 
ditions which are automatically satisfied by our non- 
perturbative approach. Therefore, the grand partition 
function Eq.(14a) can give the exact virial coefficients. 

Furthermore, the energy density and heat capacity can 
be given by the dimensionless virial coefficients in the di- 
lute Boltzmann regime. This provides a novel require- 
ment for the effective field theory. The complication 
of the temperature dependent interaction provides ad- 
ditional task to check the thermodynamic consistency 
characterized by the one-one Legendre correspondence 
relation between T and entropy S. 

According to the thermodynamical relations for s = 
S/V and e = E/V 



dP 
df 



= _ T 2JL ( ~P + Hn 
8T\ T 



and using partial derivative formulae such as 



d/j,* \ (dT\ ( d\i 



dm 2 D 
dT 



(17a) 
(17b) 

(18a) 



the integrated entropy and energy densities read 

.-id-i+i>(*=*V»' 



(19a) 



e=2fff k + + TV ( ^) ' n 2 (19b) 

J k 2m m \ m J 

In the above equations, we have employed the notation 



s ideal = - 2 / \h ^ h + (1 - fk) ln(l -/*)], (20) 
Jk 



with that for the correlation factor 



( d 2 



T>{T,n*) 



' dn\ _ d 2 n 



( dn \ 



(21) 



The first term in entropy density Eq.(19a) and first two 
terms in energy density Eq.(19b) constitute the confor- 
mal formalism of ideal Fermi gas or mean-field theory but 
with the effective chemical potential. Significantly, there 
is an additional (last) spontaneously generated term 
oc T>(T, fi*) resulting from the quantum level rearrange- 
ment. These extra terms are the explicit dynamical high 



order contributions and of non-Gaussian characteristic. 
As explicitly indicated by Eq.(21), the high order effects 
are mixed with the fluctuations through the tempera- 
ture and particle number density susceptibilities. These 
contributions are vanishing in the zero-temperature con- 
ditions. Therefore, they do not affect the previous zero- 
temperature universality property discussions[17, 21, 23]. 
Furthermore, they are canceled out by each other exactly 
in the physical Helmholtz free energy density f — e — Ts 
and consequently can be easily neglected. The present 
recipe rectifies the drawback. 

With the analytical expressions, we have indicated the 
entropy per particle versus the rescaled energy in Fig.l. 
One can see that although the high order contributions 
can be almost neglected in the low temperature regime, 
they play a significant role with the increase of temper- 
ature (energy density). In comparison to that for the 
mean-field theory without the non-Gaussian correlation 
contribution, the convex behavior of the entropy curve 
more closely approaches the experimental result or nu- 
merical simulations [6]. For comparison, we have also 
given the curve for the ideal free Fermi gas. 



E/E, 



FIG. 1: Entropy per particle versus rescaled energy at 
unitarity[6]. The solid curve is for the theoretical value in- 
cluding the spontaneously generated contribution while the 
dotted one is for the mean-field theory result. The below 
dashed curve is that for the free Fermi gas. 



B. Virial theorem and universal coefficient £(T) at 
finite temperature 

In order to justify the validity of our quasi-Gaussian 
approximation, the virial theorem for the dense and 
hot unitary Fermi gas will be strictly re-examined and 
proved. 

The non-Gaussian correction terms appear in the pres- 
sure and energy density, simultaneously. At unitarity, the 
general analytical expressions for the pressure, Eq.(14a) 
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and energy density, Eq.(19b), are reduced to 



fl /2 ( z >) f! /2 (z')f_ 1/2 (z'y 



+ 



2/1/2(^0 m /2 (z') 



= 3 - P . 
2 



(22) 



Eq.(22) explicitly demonstrates that the unitary Fermi 
gas thermodynamics obeys the virial theorem of ideal 
gas [3, 5]. What we want to emphasize is that the cru- 
cial high order non-Gaussian fluctuation terms ensure the 
virial theorem and the zero-temperature HvH theorem. 

The energy density manifests the scaling property as 
displayed by Eq.(22). The universal coefficient £(T), i.c, 
the ratio of the unitary Fermi gas energy density to that 
of the ideal ones, can be given as a function of z' 



ecn = 1 - 



Vi/2{z')h/2{z') 2f* /2 (z')fi 



C 5 /2(Z') 



(23) 



The effective fugacity z' is uniquely determined by the 
particle number density according to Eq.(8) or Eq.(16) 
for fixed temperature T. 

As shown in Fig. 2, the universal coefficient approaches 
the expected saturation value £ = 1 in the Boltzmann 
limit, ft is worthy noting that in the low temperature 
strong degenerate regime with » 1, the effective 

fugacity can be explicitly indicated as z' — e T/ / T , where 
Tf is Fermi characteristic temperature. As stressed in 
[17, 21], the zero-temperature result £(0) = | is consis- 
tent with the existing theoretical calculations [10, 12] and 
experimental measurements [6, 32, 33]. 

In Fig. 3, we have presented the finite temperature ther- 
modynamical numerical results. They are reasonably 
consistent with previous theoretical investigations [7, 8]. 
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FIG. 2: The universal coefficient £(T) versus rescaled tem- 
perature. 
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FIG. 3: a), Energy per particle versus rescaled temperature; 
b), the physical chemical potential. The solid curves are for 
the results of unitary Fermi gas. The dashed ones are those 
for the free Fermi gas. 



C. Virial coefficients in the Boltzmann regime 

To shed light on the power of the quasi-Gaussian ap- 
proximation, we will consider the virial expansion for the 
dilute unitary Fermi system. 

From the general statistical mechanics viewpoint, the 
virial equation of state is a series of pressure expanded 
in density. For the low density weak degenerate scenario 
n\ 3 < 1 , the equation of state can be expanded according 
to 



nT 



1=1 



1-1 



(24) 



The ai (T) is the border Clausius virial coefficient with g 
being the spin degenerate factor (g = 2 in this work for 
the ideal symmetric scenario). 

For the purpose of characterizing the complex corre- 
lation effects, it is worthy noting that the expansion is 
according to the effective fugacity z' = instead of 
directly according to z — e 13 ^. Meanwhile, in order to 
avoid the confusion with the conventional expansion in 
terms of z, the primes on b and c in the expansion iden- 
tities below have been explicitly indicated. In the weak 
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degenerate high temperature Boltzmann regime and by 
using such as [29, 30] 



v'2 



h/i{z') = z' - -3^ + 



z' 3 z' 4 



2 3/2 33/2 43/2 



+ ■■• 



(25) 



the pressure P and particle number density n are the 
functions of small z' and can be expanded as 



p 2 vv- 



A 3 



(26a) 
(26b) 



By comparing Eq.(24) with Eq.(26a) and Eq.(26b), the 
relations between a\, b[ and cj are derived as 



ai 



(27a) 

a 2 = 63 - c 2 , (27b) 
a 3 = 63 - 2&' 2 c' 2 + 2c 2 2 - 4, (27c) 
a 4 = 64 - 3634 + 56 2 c 2 2 - 5c 2 3 - 26 2 c 3 + 5c' 2 c 3 - c 4 , 

(27d) 

It is worthy noting that the explicit results b[ = c[ = 1 
have been used to reduce the relations. 

With the exact grand thcrmodynamical potential 
Eq.(14a) and corresponding particle number density ex- 
pression Eq.(16), one can immediately obtain the dimen- 
sionless virial coefficients at unitarity \a\ = 00 



01 = 1, 

a 2 = 0, 

as = - 
0,4 



1 4 

4 9V3 " 
15 25 



-2a 



(o) 



+ 



3 ' 
5 



32 32\/2 2^6 



_ (i + l)(i-2) (Q) 

u » — 2 » ' 



-5a 



(o) 



(28a) 
(28b) 

(28c) 
(28d) 

(28e) 



The ratios of the virial coefficients for the Unitary Fermi 
gas to their counterparts a|°^s of the ideal non- interacting 
Fermi gas are integers through a general term formula 
Eq.(28e). 



D. The second virial coefficient with non-Gaussian 
fluctuations 

From the quantum degenerate viewpoint and as in the 
classical van der Waals equation of state, the second virial 
coefficient is usually considered to be "most important" . 
Compared with a 2 °^ = ij^g of the ideal quantum gas 
with the effective quantum "attraction" (— ) for bosons 



and "repulsion" (+) for fermions, the second virial coef- 
ficient of unitary Fermi gas is found to be vanishing. 

The vanishing a 2 — is different from what was re- 
ported in the literature. In Refs.[3, 4], the simple ap- 
plication of the conventional quantum cluster expansion 
technique [29-31] gives a 2 /a 2 °^ ~ —3 (note the minus 
sign). The differences between those reported in the liter- 
ature and ours manifest on two respects: quantitatively, 
the magnitude of the second virial coefficient obtained 
in our approach is smaller than a 2 ^; qualitatively, it is 
vanishing. 

What physical reasons lead to the second virial coef- 
ficient a 2 = with the explicit difference? In applying 
the perturbative quantum cluster expansion for calcu- 
lating the second virial coefficient, the dilute condition 
n l/o 1 3 "C 1 is assumed in addition to nA 3 <SC 1. Although 
the divergent scattering amplitude should drop out in 
the final thcrmodynamical quantities, the dilute condi- 
tion n|/ | 3 < 1 is not automatically satisfied. In the 
language of resummation technique, the effects resulting 
from the infinite dynamical high orders with various irre- 
ducible expansion diagrams also contribute to the second 
virial coefficient. 

In the strongly correlated medium with n|/o| 3 ^> 1, 
especially at unitarity, the non-linear fluctuation effects 
and dynamical high order contribution have mutual influ- 
ence. The microscopic attractive dynamics and repulsive 
quantum correlations are mixed with each other. Their 
contributions for the second virial coefficient are found to 
be offset at unitarity; i.e, the competitive non-linear fluc- 
tuation effects make the expected leading order quantum 
correction to Boyle's law vanishing. 

On keeping the virial expansion up to the second or- 
der for the equation of state, the macroscopic thermo- 
dynamical property indicates an explicit Boyle tendency 
of a classical ideal gas. Therefore, the complicated non- 
Gaussian effects in the strongly correlated quantum sys- 
tem can be directly validated by measuring the first-order 
Joule- Thomson coefficient in the Boltzmann regime. 



IV. PROSPECTIVE SCALING PROPERTIES 
AND VIRIAL EXPANSION OF UNITARY FERMI 

GAS 

One the one hand, the essential divergence nature of 
the bare two-body scattering amplitude fo(k,a) = i/k 
implies that the weak coupling order-by-order "expan- 
sion" is not applicable with a real or virtual shallow 
bound state. One the other hand, the scaling property 
with the divergent scattering interaction strength implies 
that the thermodynamic expressions can be very simple; 
i.e., the analytical formulae or virial coefficients must be 
symmetrically conformal to those for the ideal noninter- 
acting Fermi gas. 

Like in the classical van der Waals equation of state, 
the lowest order of perturbative quantum cluster expan- 
sion in calculating a 2 is actually on the mean-field the- 
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ory level. In performing virial expansion at unitarity, the 
non-linear fluctuation effects should be reasonably well 
taken into account. The difficulty is reflected by the fact 
that the individual low and high order dynamical contri- 
butions cannot be separately described from the quantum 
fluctuations. 

The fluctuation effects on the thermodynamics quan- 
tities beyond the mean-field theory can be described by 
the particle number density susceptibility. Of course, the 
susceptibility itself is the final goal. With this quasi- 
Gaussian approximation method, the integrated equa- 
tion of state is uniquely determined by the particle num- 
ber susceptibility \ indirectly in terms of the auxiliary 
variable /x*, correlation factors C and V. Through /x*, 
the thermodynamical quantities can be given by a set of 
highly non-linear coupled equations. In the high temper- 
ature weak degenerate Boltzmann regime, the analyti- 
cal virial expansion can be performed by eliminating the 
auxiliary variable. 

As indicated by the coupled equations Eq.(14a) and 
Eq.(14b), the non-Gaussian fluctuation contributions 
compete with the high order dynamical effects. Their 
effects on the second virial coefficient are found to ex- 
actly cancel each other out at unitarity. Therefore, from 
the viewpoint of virial expansion to explore the quantum 
degenerate physics, one must go beyond the calculation 
of a 2 . 

A general term formula connecting the virial coeffi- 
cients of ideal Fermi gas and unitary Fermi gas is found 
for the high order virial coefficients. With the scaling 
equation of state Eq.(22), the calculated virial coefficients 
manifest the expected scaling properties. It is worthy 
noting that the negative sign in front of 04 = — 5a^ in- 
dicates a bosonization conversion tendency of the unitary 
Fermi gas at the intermediate BCS-BEC crossover point. 

V. CONCLUSION 

In summary, the three-dimensional strongly correlated 
BCS-BEC crossover thermodynamics is discussed with 
an auxiliary fluctuating mirror background shift method. 
The equations presented demonstrate a mutual influence 
of the quantum statistical correlation and dynamical ef- 
fects for different orders. 

The goal of this work is to make an analytical attempt 
doing the challenging non-linear virial series expansion. 



The obtained results manifest the scaling property of uni- 
tary Fermi gas; i.e., the virial coefficients are found to be 
proportional to those of ideal Fermi gas with a general 
formulae. 

In contrast to the investigations in the literature, the 
second virial coefficient of the unitary Fermi gas is found 
to be vanishing. To a great extent, this unpleasant result 
further puzzles us as regards reaching a deeper theoretical 
understanding of the BCS-BEC crossover physics in the 
Boltzmann regime. 

In physics, a 2 = would not be unexpected because 
the unitary Fermi system is in between the dimer boson 
and fermion phases. The non-linear fluctuation effects 
are mixed with and compete with the dynamical effects. 
Their contributions to the second virial coefficient can 
be exactly canceled out by each other. By keeping up to 
the second order virial expansion in the weakly degener- 
ate regime, the bulk properties can manifest the Boyle's 
law of classical ideal gas. The Joule-Thomson Coeffi- 
cient measurement in the unitary limit can clarify the 
difference or be used to explore the novel non-Gaussian 
correlation effects, in the affirmative. 

In the strongly interacting regime, especially at uni- 
tarity, the non-Gaussian fluctuations play an important 
role. Technically, the realistic grand thermodynami- 
cal potential f2(T, /x) or equation of state is described 
by the two coupled parametric-equations Eq.(14a) and 
Eq.(14b). Equivalent to eliminating the implicit variable, 
effective chemical potential /x*, the potential expansion 
is according to the effective fugacity z' = instead of 
directly according to fugacity z — e 13 ^ . 

The virial theorem at unitarity, the general HvH theo- 
rem at zero temperature and the fundamental thermody- 
namical consistency relations are strictly ensured in the 
formulation. This non-perturbative approach makes it 
possible to discuss the non-Gaussian fluctuation physics 
in the novel strongly correlated quantum system. 
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